Annals of the Academy of Romanian Scientists 
Series on Physics and Chemistry Sciences 
ISSN 2559 — 1061 Volume 5, Number 1/2020 38 


ALPHA-RESONANCES IN MEDIUM 
AND HEAVY NUCLEI 


I. Review of phenomenological models 


D.S. Delion*, A. Dumitrescu', and V.V. Baran? 


Abstract 


The low-lying rotational bands in light nuclei were explained a long time ago 
by a-clusters orbiting around the remaining core. These quasimolecular states 
are related to an enhancement of the elastic a-particle cross section at large 
angles, the so-called Anomalous Large Angle Scattering (ALAS) phenomenon. 
On the other hand, a-particles were evidenced by the discovery of the a-decay 
for nuclei with Z > 50. We used a phenomenological model of an a-particle 
moving in a pocket-like potential, able to explain the ALAS phenomenon, in 
order to understand the main features of the decay data between ground and 
excited states. Then we extended the Brink model, explaining giant resonances, 
and presented a possibility to detect a-like collective resonances. 
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1 Introduction 


a-clusters can be evidenced as building blocks in nuclei, due to the large a-particle 
binding energy [1, 2, 3, 4]. Let us first of all mention the pair staggering shape 
of binding energies along the so-called a-lines, defined by adding Ng = A/[2(.N + 
Z)| quartets on top of a double magic nucleus plus a given number of neutron 
pairs [5]. This behaviour is similar to the well-known even-odd staggering along 
isotope/isotone chains, induced by the pairing interaction, and lead to the concept 
of the isospin pairing between proton and neutron pairs [6]. The a-like structures 
can directly be probed in the ground states of light nuclei with a quarted structure 
like ®Be, 2C, !°O, ... 4°Ca [7]. On the other hand, clustering is favored at smaller 
nuclear densities [8], as it is the case in some excited 0* states, called Hoyle states 
with a very pronounced quartet structure [9]. 

Clustering structures were also evidenced by the a-like rotational bands in light 
nuclei, seen as resonant structures in the elastic cross section of the a-particle scat- 
tering on light nuclei. These structures, called quasimolecular states [10], are usually 
connected to the so-called Anomalous Large Angle Scattering (ALAS), seen as an 
enhancement of the ratio between the angular distribution and Rutherford distribu- 
tion at large angles. 

The a-particles were first of all evidenced by the a-decay mode for nuclei with 
Z > 50. It turns out that this process can be understood in terms of a pocket-like a- 
core interaction, similar to the quasimolecular description in light nuclei. The aim 
of this paper is to review the phenomenological approaches concerning the a-like 
resonances in a-decaying nuclei. 


2 Resonant description of the a-decay 
The dynamics of an a-daughter system emitted from a parent nucleus 

Ps Ia as (2.1) 
is described by the time-dependent Schrédinger equation 


ave 7, 22(R,t) 


oR,1) =| +ViR)| O(R,t) = ih —" (2.2) 


where yz is the reduced mass of the a-core system. The general time-dependent 
solution 


®(R,t) = U(R)e"*/" —” (2.3) 


depends on a complex energy 


B=Q-iz, (2.4) 
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where the real part denotes the Q-value of the process. The wave function describing 
a decay process should be proportional with an outgoing wave at large distances 


eikR 


where N is called scattering amplitude and « denotes the asymptotic wave number, 


defined by the usual relations 
) 
D 
eee: = 2. 
Q a p=hkt, (2.6) 


Thus, the incoming wave is missing and therefore the S-matrix of the process, pro- 
portional to the ratio between outgoing and ingoing waves, has a pole corresponding 
to this energy. By using the continuity equation, one obtains the probability flux 
being proportional to the imaginary part of the energy 


T = ho\N/? ’ v=p/p ’ (2.7) 


provided the normalisation condition in the internal region is fulfilled 


[wcrypar =1. (2.8) 


Let us mention here that the half life is defined by the standard relation 


Aln2 


For all measured a-decays the Coulomb barrier Vg = V(Rg) surounding the nuclear 
a-core interaction is large enough (Vg > 2Q), leading to the following estimate 


im 
0 < 1o-** (2.10) 


Thus, the stationary condition can be considered as a very good approximation and 
the radial wave function, satisfying the equation 


AW(R) = EV(R) , (2.11) 


decreases inside the barrier by many orders of magnitude, being very close to a 
bound state. Therefore the integration in the normalisation condition (2.8) can be 
performed only until the external tourning point. Let us mention in this context 
that the radial solution can be expanded in multipoles 


Yn (B) = 7 Ay, (A) (2.12) 
l 
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and at large distances, where the interaction becomes spherical and purely Coulom- 
bian, the equation (2.11) can be written 


a I(l+1) x 
—— “—1 = 0 2.13 
oat SS +X] fy =o, (2.13) 
in terms of the reduced radius and Coulomb parameter respectively 
p = KR 
AZ pe? 

2.14 
x ae (2.14) 


where Zp is the charge of the daughter nucleus. If the interaction potential V(R) 
is defined for all radii, the normalised wave function in the internal region with the 
outgoing boundary condition (2.5) is found for few discrete values Q; < Vg. They 
are called a-like, or Gamow resonances and at large distances are proportional to 
the outgoing solutions of Eq. (2.13) called Coulomb-Hankel waves 


filR) R400 MH} (p, x) = M[Gilp,x) + Filo, (2.15) 


written in terms of irregular (regular) Coulomb waves G (F/). The lowest eigenstate 
with k = 1 describes the a-decay process connecting ground states of parent and 
daughter nuclei with the eigenvalue equal to the experimental Q-value. This can 
be achieved by changing the depth of the a-core nuclear interaction. We will show 
later that it is possible to find in this way the first and second a-like resonances 
fulfilling the under-barrier condition Q;, < Vg. It is important to point out that 
the internal shape of the nuclear part can in principle determine the position of the 
second resonant state, which can experimentally be detected as a maximum in the 
excitation function of a-scattering cross sections close to the Coulomb barrier. We 
will analyse this possibility in a special Section. 

By using the expansion of the wave function (2.12), the total decay width be- 
comes a sum of partial components 


r=) mM? = SCT. (2.16) 
l l 


The matching condition (2.15) between each internal component and the correspond- 
ing Coulomb wave at some radius R inside the pure Coulomb barrier leads to the 
following factorisation 


Pr, = PUR)? (R), (2.17) 
in terms of the penetrability and reduced width 
20 20 
MA) = TM pt Gen 
JOM eof SEF 
he fP(R 
a(R) = EA (2.18) 


2uR 
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The theoretical decay width (2.16) is obviously different from its experimental value. 
The ratio 


Sy = bee. (2.19) 
Vth 
defines the a-decay spectroscopic factor. Let us stress here that this phenomeno- 
logical approach using the above mentioned normalisation condition of the wave 
function determined by an a-core potential cannot determine this physical quantity. 
Only a microscopic approach, describing the a-formation probability, is in principle 
able to perform this task. 


3  a-daughter interaction 


The realistic a-daughter interaction can be estimated in terms of the double folding 
between the nuclear densities [11, 12, 13], i.e. 


V(R,Q) = [een [ arapolen)palta)e(R + 2 —Tq) , (3.1) 


where R is the a-core radius, 9 are the Euler angles with respect to the symmetry 
axis of the daughter nucleus and v denotes the nucleon-nucleon force, usually taken 
as the M3Y potential [12]. This procedure was widely used to compute the potential 
between heavy ions by using for their densities a Woods-Saxon shape. In our case 
the density of the daughter nucleus is given by such a distribution, while that of the 
a-particle by a Gaussian with standard parameters [14]. 

The potential (3.1) can be divided into a spherical and a deformed component 


V(R,Q) = Vo(R) + Va(R, Q) . (3.2) 


By expanding the nuclear densities in multipoles one obtains the deformed part of 
the interaction 


Va(R, 2) = S>Vy(R)(R,Q) . (3.3) 
A>0 


Here the angular part of the wave function has the following ansatz 
IB, 2) = [Y(R) @¥A()]. (3.4) 


Due to the Pauli principle, the a-particle can exist only on the nuclear surface, 
where the nuclear density becomes less than 20% of the standard saturation value 
at equilibrium. Therefore, the local potential simulating this situation should have 
a pocket-like molecular shape centered on the nuclear surface, In Fig. 1 we plotted 
the spherical potential corresponding to the decay 7!*Po +78 Pb +a by a solid 


Alpha resonances in medium and heavy nuclei 43 


25 


V (MeV) 


R (fm) 


Figure 1: Realistic alpha-daughter double-folding plus ho repulsive potential given 
by a solid line for the system ?°°Pb + a versus radius. The repulsive harmonic 
oscillator part is plotted by dots, the parabolic approximation of the realistic nuclear 
attraction estimated by the double folding procedure by triple-dot dashes and the 
pure Coulomb potential by a dot-dashed line. The horizontal line corresponds do 
the Q-value of the emission process 7!2 Po 78 Pb + a. 


curve. The internal repulsion, given by the dotted ho potential, simulates the Pauli 
principle forbiding the existence of the a-particle inside the nucleus. The horizontal 
line denotes the energy release (Q-value). Let us mention that the three solutions of 
the equation V(R) = Q are called turning points: R,, R2, R3. Notice that the ho 
repulsion is matched to the realistic attractive part at the second turning point Ro. 


4 Universal rule for the reduced width 


Let us consider a schematic a-daughter interaction given by the ho potential, plotted 
by dots in Fig. 1, matched to a pure Coulomb repusion at Roe. By considering the 
Q-value as the first resonant eigenstate in the shifted ho well Q — V(Ro) = shw, 
together with the continuity condition at the top of the barrier Rc, one obtains the 
following relation 


Bo(Rc — Ro)? 


hwo 5 


1 
= Virag( Ro) cc 30 ) (4.1) 
where we introduced the so called fragmentation (or driving) potential 


Vyrag( Ro) = Vo(Rc) - Q . (4.2) 
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For a shifted ho well one has for the ground state 
fint(R) = Age PR F0)"/? (4.3) 
By using Eq. (4.1) one obtains the following relation 


2 h 2 loge? 2 ‘ 
logio ¥ (Ro) = logo IpRe tint(Re) = Fw Virag( Rc) + logig QeuRo , 


0 
(4.4) 


called the universal rule for the reduced width [15]. Let us mention that it does 
not depend upon the equilibrium radius Ro and the ho frequency of the internal 
potential hwo can be considered as a free parameter, together with the matching 


radius Ro. 
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Figure 2: Logarithm of the reduced a-decay width for even-even emitters versus the 
fragmentation potential Verag = Vo(Rc) — Q, for Ro = 1.2(A73 + 41/3). By open 
symbols are given emitters above 7°°Pb and by dark symbols all other emitters. 


In Fig. 2 we plotted the logarithm of the reduced width estimated at the geo- 


metrical touching radius 


Ro = 1.2(Ag? + 41/3) (4.5) 
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by considering in panel (b) even-even emitters with 82 < Z < 102, N > 126 and in 
panel (a) all other cases. By open symbols are given emitters above 7°°Pb and by 
dark symbols all other emitters. One clearly sees that the slope is indeed negative, 
as predicted by the universal law (4.4). This is clear evidence of the a-clustering on 
the nuclear surface, seen as an a-like resonace in a pocket-like potential. 


logiol (MeV) 
4 
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Figure 3: Logarithm of the a-decay width for even-even emitters versus the frag- 
mentation potential log P + aVfrag + b. The symbols are the same as in Fig. 2. 


In Fig. 3 we plotted the logarithm of the decay width versus log,y P+aV frag +8, 
where we considered different fitting parameters for the two panels in Fig. 2, namely 

(a): ay = —0.353, by = +0.028, hwo = 17.014 MeV; 

(b): ag = —0.139, bo = —3.013, hwy = 43.055 MeV . 

This plot expains the well known universal Viola-Seaborg dependence [16] versus 
the variable V = (q.Zp + c2)Q71/? + c3Zp + c4, due to the fact that logig P ~ 
(c1Zp + @jQ-/? and Verag ~ Zp. 

Let us finally mention that this universal law for the reduced width is also valid 
for proton and heavy cluster emission [15]. 
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5 Coupled channels approach for a-like resonances 
For the emission process connecting some states 
P(Jp) > D(J) + a(Z), (5.1) 


the wave function describing the a-daughter dynamics can be expanded as follows 


N 
0), (R) => 5 (0, &), (5.2) 


where c = (J, L, Jp) is the channel index and the core-angular harmonics, depending 
on the daughter coordinates 2 and a-particle angles R, are defined by [17] 


y(2,R) = [WP Q evi). (5.3) 
JpMp 
For a-decays from the ground states (gs) of even-even nuclei with Jp = 07, one has 
c= J = L=even 


y1(0, R) = [¥}P)(Q) @ ¥1(8)] (5.4) 


In the case of rotational nuclei we have vu) (Q) = Yzm(Q), and the final state 


corresponds to a “quasimolecular state” in which the daughter nucleus and a-particle 
rotate in opposite directions with total momentum Jp = 0. For vibrational nuclei, 
yw) is a vibrational wave function and the picture corresponds to a “quasimolecule” 
where the two partners vibrate in opposite directions. 

For transitions from odd-mass nuclei, Jp and J have half-integer values and the 
a-particles are emitted with angular momenta |Jp — J| < L < Jp+J. The resonant 
states have positive parity for L=even and negative parity for L=odd. The daughter 
wave function has a core-quasiparticle ansatz 


wy (Q, r) LX $7(Q )@ oj(r Visas» (5.5) 


where the coefficients are found by diagonalizing a core-quasiparticle interaction 
reproducing the energy spectrum. Notice that a rotational band with J > 7 is 
defined by considering only a given quasiparticle orbital 7 in the above superposition. 
Let us mention that the Coherent State Model (CSM) description describes the 
two limits in a unified way by changing a “deformation parameter” d, which enters 
the coherent superposition of quadrupole bosons exp(db! xe) [LS 19]. 
The a-daughter dynamics is described by the stationary Schrédinger equation 


AW 7,up(Q, R) = QU z,u,(, R). (5.6) 
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An a-decaying state is identified with a narrow resonant solution that contains only 
outgoing components, called Gamow resonances. The Hamiltonian 


h2 
H= a Ve + Hp(®) + Vol, B) + Va(Q, R) (5.7) 


contains the kinetic term, a term describing the dynamics of the nuclear core Hp(Q) 


D D 
Hwy, (2) = BV}, (O), (5.8) 
and the a-core potential, which we split into spherical and deformed parts. By 
using the orthonormality of the core-angular harmonics in the superposition (5.2), 
one obtains the system of differential equations for radial components [17] 


N 


= Ad he2) 6S lag (5.9) 


c=! 


dp(R) 
dpz 


where the matrix 


(c) (c’) 
Le(Le +1) | Vo(O, RB) (V7; 1Va(Q, RV,” 
Ace! = -—$——$ So 1 cc! a . 
(8) | pa 7 Qa _ E ° 7 Qa a E, ° yy 
is given in terms of the channel reduced radius 
2 Ons E, 
Po = KR, Ke= ae (5.11) 


Let us mention that at large distances the potential becomes spherical (Vz — 0) 
and purely Coulombian. The external solutions in each channel c can be found by 
superposition 


N 
perO(R) = SH RIN. (5.12) 
a=1 


Here we have involved the columns of the fundamental matrix of Coulomb solutions 
with outgoing Coulomb-Hankel asymptotics 


HD (R) Gea(R) + iFea(R) ———> SeaH\" (R) 
+00 ec 


= Og (Gin) +i {RA )). (5:13) 
where c denotes the channel and a the resonant eigenvalue index. Inside the Coulomb 


barrier, it coincides with the real matrix of irregular Coulomb solutions HS? (R)| = 
Gea(R). The unknown coefficients NV, are called scattering amplitudes. 
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From the continuity equation one obtains the total decay width as a sum of 
partial widths [17] 


P= DoPe= } 7 hve lim | fo(R)|? 


S > hive Nel”, (5.14) 


depending on scattering amplitudes and asymptotic velocities uv. = K-/y. By invert- 
ing Eq. (5.12) and using (7.21), the external components of the wave function at 
some point R inside the a-daughter potential can be expressed 


N 
T 
(cat) _ (H(p),/—™ cx 
per) (R) DH (Ry pe = Lew N, (5.15) 


in terms of the partial decay widths [, and asymptotic channel velocities. 

Let us mention that for the spherical one-channel limit with L.=0, by considering 
the matching condition y@) (R) = W*(R), Eq. (5.15) becomes the well known 
relation 


v6 (RB) ]” 


Goo(R) oe 


To = hvo | 


6 Systematics of a-clusters on the nuclear surface 


Microscopic calculations within the mean field plus pairing approach evidenced that 
the a-formation amplitude has a shifted gaussian shape centered on the nuclear 
surface [17], mainly concentrated in the monopole component. This conclusion is 
also supported by a systematic analysis of particle and heavy cluster emission pro- 
cesses [15]. The coupled channels analysis confirmed a similar behavior of resonant 
solutions [19] 


; 1 

yh")(X) = Aoyf yf er *7? (6.1) 
T 

where X = R— Ro is the shifted radial coordinate and N the number of channels. 

Here, Ap denotes the cluster amplitude and (o the ho parameter of the monopole 

component 


1 MoWo 
= —— = = hw . 
Bo be i f 0 
2 
f = [5 ~0.096 (MeV-! fm-), (6.2) 


(he)? 
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which can be written in terms of the ho length parameter bo. We considered that 
for heavy nuclei one has uw * ma. This wave function is generated by a shifted ho 
a-daughter local potential 


in 1 
Voi") (X) = up + 5 liw0/5oX°. (6.3) 


We matched this potential to a realistic external interaction, given by the dou- 
ble folding a-nucleus potential in terms of the M3Y nucleon-nucleon plus Coulomb 
proton-proton interaction [19]. The parameters were determined by a-scattering 
experiments. For the a-particle, we assumed a gaussian distribution with bg=1.19 
fm [11, 12]. 

We match the internal shifted ho potential to the monopole component of the 
external double folding potential Vo(R). The same remains valid for the wave func- 
tion. The unknown coefficients vg, G9 and Ro, are determined by using the following 
matching relations 


1 
vo + 55 AOXim = VW(Rm), Xm=Rm— Ro 

1 

70% = Vo(Rm) 
in (Xm) = —BoXm = In’ Ye (Rm). (6.4) 


The Q-value is given by experimental data and therefore one obtains the ho param- 
eter and X,, from the last two equations 


In’ be") (Rm) 
— wn VE) (Rn) | FV (Rm) 
cc “Wace ies ae (6.5) 


We fix the matching radius at the second internal turning point R,, = Re, thus 
ensuring the existence of a narrow resonance corresponding to the first eigenvalue 
in the ho pocket with Q — vp ~ sliwo. One determines (9 and X,,, and therefore 
the equilibrium radius Rp. The value of vp is given by the first equation (6.4). The 
internal channel wave functions are given by gaussians with slightly different ho 


parameters 
: 1 Be = 
per) (X) = A.\| —1/—e BoX?/2. (6.6) 


The ho parameters and amplitudes are determined by the standard matching con- 
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ditions 


_ mo 1 (ext) 
Be —_ Xu In We (ae) 


Ac = VE (Rin), Ne (6.7) 


The a-particle spectroscopic factor is given by 


Sq = / |U(R)PdR~S > AZ = SoS. (6.8) 
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Figure 4: Vibrational resonant states (solid lines) and rotational bands built on them 
(dotted lines) in the a —° Ce molecular potential. The lower solid line corresponds 
to the Q-value of the process. The dot-dashed line is 10 times the ratio p(R)/po and 
the dashed line shows a factor 20 multiplying the gs a-particle probability versus 


radius [20]. 


We first analysed the resonant eigenstates in the pocket-like a-core potential. In 
Fig. 4 we plotted a typical spectrum of resonant states for the transition 44Nd —'40 
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Ce-+a [21]. The first resonant zero node state has the eigenvalue Eo corresponding 
to the Q-value, drawn by the lower solid line. Lower solid and dotted lines are 
eigenstates in the parent nucleus with Jp = 07, described by coupling the low-lying 
eigenstates of the daughter nucleus with a-particle states [ZL ® L]p in Eq. (5.4). 
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Figure 5: Monopole ho length parameter versus neutron number of the daughter 
nucleus (a) and versus the ratio between the energy difference of vibrational levels 
and fragmentation potential (b) for even-even emitters [20]. The symbols denote 
different regions of the nuclear chart, divided by magic numbers. The regression 
line in panel (b) fits the data, excepting the first region (dark squares). 


We considered experimental energies of the daughter nucleus in the coupled sys- 
tem of equations (5.9), (5.10). A systematic analysis of their B(E2) values is given 
in Ref. [19]. In Ref. [22], we have shown that for 7!?Po the low-lying 0+ a-like reso- 
nances (except gs) have small spectroscopic factors and they can be experimentally 
detected only above 6 MeV. We expect a similar behavior for low lying monopole 
resonances, plotted by dotted lines in Fig. 4. Thus, they cannot be detected ex- 
perimentally and the 0* states observed around 2 MeV in '44Nd have a different 
nature. 

The upper solid line denotes the next eigenvalue E, of the first one-node excited 
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Figure 6: Monopole ho length parameter (a) and the ratio between the energy 
difference of vibrational levels and fragmentation potential versus neutron number 
of the daughter nucleus (b) for odd-mass emitters. 


vibrational state and the dotted lines above it belong to the excited monopole band 
with Jp = 0*. The decay width [9 in Fig. 4 is given by the barrier penetration 
Pren = 5 10-4! keV, multiplied with the spectroscopic factor Sq = 0.125. It repro- 
duces the input experimental value Pez, = 6.3 10-4? keV. Moreover, our numerical 
estimate gives a ratio between errors 61 /dF,=0.16. 

The energy difference AEF = E, — Ep determines the ho parameter of the molec- 
ular potential. We plotted the wave function multiplied by 20 with a dashed line, 
and with a dot-dashed line the ratio between the nuclear and equilibrium density 
multiplied by 10 

p(R) 1 


fT Pee RS a (6.9) 


by considering a standard diffusivity a = 0.5 fm and nuclear radius R,, = 1.2 A?/®. 
These values are close to the systematics of Ref. [23] based on electron scattering 
data. Notice that the Mott phase transition point between the nucleonic and a- 
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clustering phases corresponds to 10% of the equilibrium density in infinite nuclear 
matter [8] and 20% for 7!*Po, as seen in Fig. 5 of Ref. [24]. Therefore our calculation 
confirms that the maximal value of the wave function corresponds to a radius larger 
than the phase transition point. 


As we mentioned before, the rotational quasimolecular band is connected to the 
ALAS phenomenon. In Ref. [25] it was pointed out that the experimental evidence 
of this effect was found for a lighter configuration with Z = 40, namely a+%°Zr. 
Anyway, one sees from Refs. [26, 27] that for the configuration analysed above, 
a+Ce with Z = 58, the ALAS phenomenon disappears close to the Coulomb 
barrier Vo ~ 16 MeV, due to the fact that exchage effects diminish for Z > 50. 


In Fig. 5 (a), we plotted the ho length parameter bo of the gs channel (L=0) as a 
function of the neutron number in the daughter nucleus for 162 even-even emitters. 
Notice that bo € [0.6, 0.8] fm, i.e. the ratio to the standard a-particle length param- 
eter by = 1.19 fm belongs to the interval bo/ba € [0.5,0.7]. In panel (b) is given the 
length parameter versus the ratio between the energy difference of vibrational levels 
AE = E, — Eo and the fragmentation potential Vfrag = Vo(Rmaxr) — Q. Notice that 
Virrag € [10,17] MeV. Except for the first region, Z ~ N ~ 50 of “superallowed” 
a-transitions [28], these quantities are linearly correlated. For those nuclei with 
bo > 0.75 fm, the rotational band built on the first excited state lies close to the 
Coulomb barrier AE'/Vyrag ~ 1 and therefore it can be detected as a structure of 
peaks in the energy dependent cross section. The regression line in Fig. 5 (b), fiting 
the data excepting the first region (dark squares), has a mean error (dbo) = 0.03 fm. 


The fine structure of odd-mass emitters was measured for Z € [87,100]. In 
Fig. 6 we show plots similar to Fig. 5 for odd-mass emitters. Here, one obtains 
bo € [0.6, 1] fm. AE/Verag € [0.4, 2]. Additionally, the L = 1 angular momentum of 
the emitted a-particle is allowed. Thus, the lower and upper vibrational eigenstates 
in the pocket-like potential of Fig. 4 can be connected by a dipole operator D ~ 
X = R— Ro of the a-particle center of mass. Therefore, the second vibrational state 
can be detected in the excitation function of an incident y beam. Such beams are 
produced at operating y-beam facilities, and will be provided in the future at the 
ELLNP facility [29]. 

In Fig. 7 we plotted the spectroscopic factor versus the neutron number of the 
daughter nucleus. Let us stress on large values above the magic numbers Z = 50, 82 
and N = 82, 126, where a-molecules are born with significant probabilities [19] and 
the best candidates for experimental detection have AE /Verag ~ 1. 


7 Resonances in the a-decay fine structure 


Transitions to excited states are a sensitive tool to probe nuclear structure details 
and to evidence a-clustering on the nuclear surface. 
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Figure 7: Spectroscopic factor versus neutron number of the daughter nucleus for 
even-even emitters. The symbols denote different regions of the nuclear chart divided 
by magic numbers [20]. 


7.1 Coherent State Model 


In Ref. [19] a systematic description of electromagnetic and a-transitions of even- 
even nuclei was performed in terms of the Coherent State Model (CSM) proposed 
bty Ref. [30]. This is a tool describing in a unified way the spectra of vibrational, 
transitional and rotational nuclei. It treats surface vibrations of deformed nuclei by 
using an exponential superposition of boson operators. The model was later devel- 
oped for the description of low-lying as well as high spin states in nuclei, including 
isospin degrees of freedom (for a review, see Ref. [31}). 

The wave function of an axially deformed even-even nucleus in the intrinsic 
system of coordinates is given by a coherent superposition of quadrupole boson 
operators bo,, with 4 = 0 acting on the vacuum state 


by) = eH) 1g) , (7.1) 


in terms of the deformation parameter proportional to the static quadrupole defor- 
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mation 
d= KB . (7.2) 


Physical states which define the ground band are obtained by projecting out the 
angular momentum 


eM) = NG Parola) (7.3) 
in terms of the projection operator 


a 2J4+1 a 
Bin = gar fw Dirx(w) Rw), (7.4) 


where Dj, ,-(w) is a Wigner function and R(w) is a rotation operator, parametrized 
by the Euler angles w. The simplest estimate of the ground band energy spectrum 
is given by 


By(d) = Ar [oP Mp9) — (PN ley?)] (7.5) 


where N is the operator for the number of bosons. Notice that for small values of 
d, the energy spectrum has a vibrational character Ey ~ A,J, while for large values 
it has a rotational shape Ey ~ A, J(J +1) [30]. 


7.2 Coupled channels description of a-transitions 
in even-even nuclei 


We apply here the coupled channels framework previously discussed, by considering 
that the daughter wave function is described by the CSM. The wave function of the 
a-daughter system has the total spin of the initial ground state (i.e. zero) 


> I) 2 (62,9) 
J 


le? e ¥3(9)] (7.6) 


WV (be, R) 


Z (bo, Q) 


Here, of) is the J-th eigenstate (7.3) of the CSM Hamiltonian Hp(b2) in terms 
of the quadrupole boson bg describing the dynamics of the daughter nucleus and 
R = (R,Q) denotes the distance between the centers of the two fragments. The 
a-daughter dynamics is described by the stationary Schrddinger equation 


HW (b2,R) = EWV(b2,R) . (7.7) 
The Hamiltonian describing a-decay is written 


h2 
A= rad: + Hp(b2) + V(b2, R) (7.8) 
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in terms of the reduced mass of the dinuclear system p. We estimate the interaction 
between nuclei as a sum of two terms 


V (b2,R) = Vo(R) + Va(b2,R) , (7.9) 


where the monopole part of the interaction is given by the same ansatz as in Ref. 
[32], i.e 


Vo(R) = VaV 0(R) » R>Rn 
= e(R-Rmin)—u, R<Rm- (7.10) 


Here, Vo is the nuclear plus Coulomb interaction, estimated by using the double 
folding procedure within the M8Y particle-particle interaction with Reid soft core 
parametrisation [11, 12, 13]. The value vg = 1 corresponds to a “pure” a-cluster 
model. By considering vg < 1 one assumes an a-cluster probability less than unity, 
necessary to reproduce the experimental half life. The second line is the repulsive 
core simulating the Pauli effect and fixing the energy of the first resonant state to 
the experimental Q-value. 

We applied the procedure of Ref. [32] to determine the matching radius R,,, and 
the coordinate R,,;n,, corresponding to the minimal value, by using the equality be- 
tween the external attractive and internal repulsion, together with their derivatives 
(see Eqs. (32) of this reference). Thus, the above interaction is continuous and it 
depends upon only one independent parameter, due to the fact that the repulsive 
strength c is inversely proportional with respect to the potential depth vo [32]. 

The \ = 2 term is given by the quadrupole-quadrupole (QQ) interaction 


dVo(R) 
dR 


Vile. R) = -Qih- Aa) 2[Q2 @ Yo(Q)Ip - (7.11) 


By using the orthonormality of angular functions entering the superposition (7.6) 
one obtains in a standard way the coupled system of differential equations for radial 
components [17] 


FIR) =e )fr(R) , (7415) 
p5 


where the coupling matrix is given by 


J(J +1) Vo(R) 
, = ———_—_—___- _- 1 / Z Zs £ 
Ajy(R) a +E Jay + By g|Vo(b2,R)|Zs") , 
(7.13) 
in terms of the reduced radius 
2u(h -— Bb 
py = KR, KZ = dear (7.14) 


h2 
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The matrix element of the particle-core coupling entering Eq. (7.13) is given in Ref. 
[30] and it is proportional to the reduced matrix element defining electromagnetic 
transitions, but with a different anaharmonic parameter 


_ —\d¥o(R) 1, @ (9) 
(Z3|Vo(be,R|Zy) = —Co(R—- Rmin) IR ATR es lQ2|ley’) (Ya||¥al| Yor) 
AV(R) d J ie 
= -—O(d)(R—-Rmin —— = (J; 20|.J'0 
(a) ( ) Tae Fr 0 20"0) 


Te. dN fas) 
ING INO 


where we defined the effective a-daughter coupling strength 


C(d) =Co ( = i Fon) (7.16) 


In order to solve the system of quations (7.12) let us first define the internal and 
external fundamental solutions satisfying the following boundary conditions: 


Rj1(R) ROH sare s , 


[o,e) 


HY (R) =Gy(R)+tFy(R) “29 by,HSP (eR) = 551 [Ga(esR) + iF s(KsR)] 
(7.17) 


Here, Ro is a radius inside the internal repulsive potential and ¢, are arbitrary small 
numbers. The index J labels the component and IJ solution, Gy(K,;R), Fy(K7R) 
are the irregular and regular spherical Coulomb wave functions depending on the 
momentum «7 in the channel J. 

Each component of the solution is built as a superposition of N independent 
fundamental solutions. We impose the matching conditions at some radius R, inside 
the barrier 


fx(fi) = S- Rar Ri) Mr = SHG? (Ri) Nr 
T T 


dfs(Ra) dR s1(R1) dH Sy (Ra) 
TR > dR p= I (7.18) 
where Ny are called scattering amplitudes. The conditions (7.18) give the following 
secular equation 


R(Ri) HO (Ri) 
dR(Ri1)  dH)(R1) 
“dk dR 
R(Ri) G(R) 
dR(Ri) dG(R1) 
“dr © dR 


=0. (7.19) 


aw) 
~ 
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The first condition is fulfilled for complex energies, determining the resonant states. 
They practically coincide with the real scattering resonant states, due to the fact that 
the imaginary parts of energies are much smaller than the corresponding real parts, 
corresponding to vanishing regular Coulomb functions F’; inside the barrier. The 
roots of the equation (7.19) do not depend upon the matching radius R;, because 
both internal and external solutions satisfy the same Schrédinger equation. The 
unknown coefficients M;,.N; are determined from the normalisation of the wave 
function in the internal region 


R3 
> [iter =1, (7.20) 
ys 0 


where R3 is the external turning point. All known half-lives in @ emission are much 
larger that the characteristic nuclear time Tyin © 10~° s >> Ty & 107-2? s. Thus, 
any a-decaying state practically behaves like a bound state, having an exponential 
decrease versus radius inside the barrier. 

From the continuity equation one obtains the total decay width as a sum of 
partial widths 


T — Doty = Do fies Jim | fs(R)/? 
J J 


S > hus |Nol? ; (7.21) 
J 


where vz is the center of mass velocity at infinity in the channel J, i.e. 


_ hy 


: (7.22) 


Us 


We used the CSM deformation parameter d, given by analyzing low-lying energy 
levels, in order to determine the a-decay fine structure, defined by the logarithm of 
the ratio between decay widths to ground and J™ states 


i 
Ty = logig = ; (7.23) 


where partial widths are given by Eq. (7.21). We call this quantity the intensity of 
the a-decay to the J-th state. In our calculations we considered vg = 1 and fixed the 
value of the repulsive strength c = 50 MeV, as in Ref. [33], thus changing the only 
remaining free parameter vg in order to reproduce the Q-value for each a-decay. 
By using the effective a-daughter coupling strength C(d), we reproduced the 
available experimental quadrupole intensities Jz. In Fig. 8 (a) we plotted the values 
of this strength as a function of the deformation d. Thus, we indeed obtained the 
linear dependence of the a-daughter QQ coupling strength with a negative slope, 
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Figure 8: (a) a-daughter QQ coupling strength C defined by Eq. (7.16) versus the 
deformation parameter d. (b) Same as in (a), but for the mass number A [18}. 
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Figure 9: (a) a-daughter QQ coupling strength C defined by Eq. (7.16) versus 
the a-particle reduced width. (b) Same as in (a), but versus N — Nmagic, where 
Nmmagic = 126 for N > 126, Nmagic = 82 for N < 126 [18]. 


predicted by Eq. (7.16). Notice that for the region N < 126 the values belong a 
narrow interval of small values C’ € [0.05, 0.1]. 


60 D.S. Delion, A. Dumitrescu and V.V. Baran 


It turns out that the QQ strength is also related to the mass number A. In 
Fig. 8 (b) we give the dependence of this strength versus A. The a-daughter QQ 
coupling is stronger for nuclei above ?°°Pb, where C * 0.5, i.e. in the region where 
the a-clustering is larger [15]. It decreases by one order of magnitude around the 
mass region A = 240. By crossing the shell closure to the region N < 126 it jumps 
by one order of magnitude, up to the value C' = 0.05. Let us stress here that for 
N > 126 this behaviour is similar to the dependence of the reduced width (called 
a-particle probability) for transitions connecting ground states with J=0. Indeed, it 
turns out that the a-daughter strength is proportional to the a-particle probability, 
with different slopes for N > 126 (dark circles) and N < 126 (open symbols), as we 
can see in Fig. 9 (a). Notice that the dependence versus the number of neutrons 
above the closed shell Nmagic is common for both regions, as we can see in Fig. 9 
(b). This proves the a-like resonant quasimolecular picture of an a-cluster moving 
around a doubly magic core. Several a-clusters are rapidly “disolved” due to the 
Pauli principle. 
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Figure 10: (a) Predicted intensity 4, by using the coupling strength C(d) reproduc- 
ing the corresponding value Jj, versus the index n in the first column of the Table 
I. (b) Same as in (a), but for Jg [18]. 


We used these values of the a-daughter strength C(d) reproducing J values 
in order to predict the intensities J, and Ig. They are plotted in Fig. 10 by dark 
symbols versus the index n. These values reproduce the available experimental data, 
plotted by open circles, with a reasonable accuracy. This proves that this relatively 
simple model is able to simultaneously describe all available experimental a-decay 
intensities to excited states within a range of four orders of magnitude. 
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8 a-like collective resonances 


We will investigate the possibility of a new collective motion performed by a-clusters 
with respect to the remaining core [34]. We will analyze the collective modes gen- 
erated by the relative motion of protons, neutrons and a-clusters, by generalizing 
the procedure proposed by D. Brink to describe giant resonance in Ref. [35] and 
used in Ref. [36] for the pygmy resonance. Let us consider four subsystems: Z, 
core protons, N. core neutron, N, excess neutrons with the same mass m and Q 
a-clusters with effective mass ma, moving in a ho mean field 


i=1 i=1 a=1 i=1 


The frequency wo is obtained from the requirement to reproduce the mean square 
radius of a sphere of nuclear radius, its value being hw) = 41 A~!/2MeV. This 
Hamiltonian may be arranged in the following form 


2 
Phot Pe Pe 
2nZ. 2N. ZN. 2wNaQ 


2 2 2 2 
Mw, Mw, Mu, MaW 
- ot °Z.Re + a CN =e ON.R? + . °QR2 


+ internal part . (8.25) 


Let us define the core, core+excess neutrons, effective a-cluster and total mass 
numbers 


Ma 


Asa Geena, ApsAeN.,. Apes —". ASApPAs. 320) 
m 


The collective part of the Hamiltonian becomes 


1 /P2 p2 p2 p? mu? 
a (F a tae “) sz (ZeRp + NeRn + NeRe + AoRa) (8.27) 
c c Ee (ey 


Let us separate the center of mass (cm) motion of the nucleus from the other modes 
of oscillation. The cm position is 


ZR, + NeRy + NeRe + AgRe 


8.28 
Ze+Ne+Ne+ Aa ( ) 


Rom = 


while the coordinate describing the oscillation of the core protons against the core 
neutrons is 


Xc=R,—-R,. (8.29) 
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The motion of the core against the excess neutrons is given by 


Z-Ry + N-R 
Y = Remyore — R. =—— —— 


= ee Ry 8.30 
OS te Ne, Z eon) 


We now introduce the coordinate which gives the motion of the system core+excess 
neutrons against the clusters 


— ZRptNeRnt+NeRe _ 


Yo _ Rem, core+exc = Ra = Zn a Ne she Ne R, . (8.31) 
The Hamiltonian acquires the diagonal form 
p2 Mywe 
A= — oxX% 8.32 
Ss (ses tex) | (2.32) 


K=CM,C,y,a 


where the effective masses corresponding to each mode are given by 


NZ, NA AqAo 
Mom =™mA, Mas M, =m Ay” M,=m 7 . (8.33) 
Let us now introduce the dipole operator 
Ze Q 
D = e) (tpi-Rom)+€o > (Tai — Rom) 
i=l i=1 
= eZ-(Rp _ Rom) + Ca Q(Ra _ Rem) i (8.34) 


We use the notations for the ratios of the effective charge and mass of the cluster 
to those of the nucleon € = e,/e, 4 = M_a/m. By adjusting the values of these 
parameters, the effects of clusters of any size may be evaluated. However, the 
spectroscopic factors for heavy clusters are much smaller than in the a-cluster case, 
making such oscillations unobservable [37]. 

The dipole operator is divided into terms belonging to (1) the core-proton against 
core-neutron oscillation, (2) the core against extra-neutron oscillation and (3) the 
oscillation of the alpha particles against all other constituents, i.e. D=Do+Dy,-+ 
Da, where 


ZN, 
Do = g “Xo es dcXc > 
Ae 
ZN, 
D, = ras — dyY 5 
Ag € 
Dy = 2 (Z.==Ap Ye Saaves (8.35) 
A LL 


The residual dipole-dipole interaction to the Hamiltonian of Eq. (8.32) is neces- 
sary in order to reproduce the experimental values of the GDR frequency. Previous 
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studies investigated the coupling of the X¢ and Y modes by the dipole-dipole inter- 
action in schematic models based on the on the Tamm-Dancoff and random phase 
approximations [38] and within the ho shell model [39]. In the present work, we 
generalize the latter by including the mode corresponding to the a type resonance 
Yo 


il 1 
Age = 5 xkKL Dk - Dr = 5 > xKL dkd, XK-Xp. (8.36) 
K,L=Cy,y,a K,L=Cyy,a 


The Hamiltonian AH; of our model is thus obtained by taking into account Eqs. 
(8.32) and (8.36), Hiop = H+ Hinz. The normal modes are obtained as superpositions 
of the individual modes C,y,a, after the diagonalization procedure. Each normal 
K mode has a dominant contribution from the original Cy, or @ modes, and the 
total Hamiltonian becomes Hig = Ue K=1,23 Ar. 

The total dipole absorption cross section of a given normal mode K, giving the 
fraction of the energy weighted sum rule (EWSR) exhausted by the mode, may be 
estimated by using the Thomas-Reiche-Kuhn sum rule (see standard textbooks, e.g. 
Appendix B of [40], Section 5.6 of [41]) 


Ane? 


Are? 1 
= = 
- he 


DE — Eo) |(f|Dx\0)? = ==> 


(0| (Dx, |Hx,Dx]]|0) (8.37) 


where Dx is the dipole operator corresponding to the mode K. In our schematic 
approach we consider only the dipole-dipole residual interaction, but it is important 
to mention that other components of the nucleon-nucleon interaction, especially the 
tensor force, are relevant for the full understanding of the dipole strength function 
and the energy-weighted sum rule. 

We simulated the nuclear mean field by a schematic ho well neglecting the spin- 
orbit interaction. In cluster models, where composite objects are formed from nu- 
cleons lying on the mean field levels, the Pauli principle may be taken into account 
by restricting the values of the radial and angular momentum quantum numbers of 
the cluster component. The antisymmetrisation effects are reproduced if the num- 
ber of cluster excitation quanta G = 2N + L is restricted to values greater than 
the sum over the single particle quantum number of the cluster nucleons that are 
excluded from below the Fermi energy, }>.(2n¢ + /<). This approach, known as the 
Wildermuth rule, is widely used to describe alpha and heavy cluster formation in 
a potential well, by using an effective value of the quantum number G in order to 
satisfy the Bohr-Sommerfeld rule entering the quantization condition [42]. In our 
case, the nucleons and clusters move in the same ho potential. The dipolar operator 
promotes the clusters on a higher level, changing the angular momentum/parity by 
one unit. Thus, all nucleons forming the clusters are automatically promoted above 
the Fermi sea, satisfying the Wildermuth rule. In a microscopic picture, they are 
localizated on the nuclear surface, and the cluster state may be thought of as an 
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additional Gaussian attractive component on the surface region of the mean field 
(28, 43]. The overlap of the cluster state with the internal core functions is negligible, 
as a manifestation of the Pauli principle. In our model, the only restriction is con- 
nected to the strength of the dipole-dipole interaction, which should have a definite 
value. Obviously, for a realistic Woods-Saxon mean field the excited resonance is 
found in the continuum and practically any value of the dipole strength is allowed. 

In order to analyze the interplay between the usual giant dipole resonance (GDR) 
and a-like resonance we will first consider the symmetric system with N = Z, where 
the pygmy resonance induced by neutrons in excess disappears. The fraction of the 
EWSR exhausted by the a-like mode may be estimated from Eq. (8.35), as follows 


(8.38) 


Figure 11: a-mode EWSR fraction vs effective cluster mass and number of clusters 
[34]. 


The alpha clusters do not exist with a probability of unity inside the nucleus. 
This fact is quantitatively described by the cluster spectroscopic factor, estimated 
as the spatial integral of the a-particle formation amplitude squared. It has the 
maximal value S, < 0.01 for lightest a-emitters above 1°°Sn [28]. For instance, a 
spectroscopic factor Sq = 0.01 corresponds to the case where for 1/100 of the time 
one has a real alpha-particle with four nucleonic masses and 99/100 of the time one 
has free nucleons. We simulate this probability by an effective alpha-particle mass. 
Thus, we define the effective mass by the probabilistic average 


Ma = PaMa,o0 + (1 = Dex tT ; (8.39) 


where m, is the nucleon mass, ™q,9 = 4M is the bare a-cluster mass, and po = Sq 
is the probability of formation of the a-cluster. 
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Therefore one has 4 < 1.03 and the smallest value of the a-mode EWSR fraction 
is fa ~ 10% above the the doubly magic nucleus !Sn. An interesting conclusion is 
that for light systems, where yu ~ 4, this fraction has lower values. This behaviour 
is compensated by the increase of a-cluster number, as can be seen in Fig. 11, 
where we plotted fa versus the number of clusters for various effective masses. This 
quantity increases versus the number of clusters but decreases by increasing the 
spectroscopic factor (effective mass). 


Let us consider the influence of the pygmy resonance induced by the excess 
neutrons. We will work within the zero inter-mode coupling case. By allowing the 
presence of excess neutrons, the a-fraction of the EWSR becomes 


2 
| ea 


The EWSR exhausted by the pygmy mode may be computed in a similar manner. 
The ratio of a and pygmy modes’ exhausted EWSR is given by: 


2 
2€ A. 2 /N. 
LS pl Al 
( *) Ne pV Ac ea 


We plot the dependence of this ratio as a function of the effective a-mass and the 
excess neutrons in Figs. 12 and 13. 


fa _ wQ 


fy A 


Figure 12: Ratio between a and pygmy EWSR fractions vs excess neutrons and 
effective a mass (Q=1). Only values larger than unity are shown [34]. 
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Figure 13: Ratio between a@ and pygmy EWSR fractions vs excess neutrons and 
effective a mass (Q=3). Only values larger than unity are shown [34]. 


The frequencies of the GDR, pygmy and a modes may be estimated as follows 


2 2 Xe Ac 
SEDER: = 0 a 4 
Be ne ae Ne. 
y Om 4 Ae + Ne? 
A QA 2€ NEY" 
2 2 C c e 
= 1-—-—{,1+— 8.42 
Wa “0 T Xa‘gn A(A. + Ne) | ( + =) ee?) 


In particular, in this case of decoupled modes, one obtains 


2 2 
Wo ~ 9 _ Xo ta 
we — we Xy Fy 


= (8.43) 


The a@ and pygmy frequencies are very close to each other, around 8 MeV above 
100Sn. In Fig. 12 we plotted the ratio between a and pygmy EWSR fractions versus 
the excess of neutrons and effective a mass for the number of quartets Q=1. As 
expected, this ratio decreases by increasing the number of excess neutrons. As we 
can see, the unexpected decreasing behaviour versus the increase of the effective 
mass (i.e. spectroscopic factor) is similar to the previous N=Z case. In Fig. 13 we 
plotted the same dependencies, but for Q=3 a-clusters. One sees that the results 
are very close to the previous case. 

Let us mention that this model has no theoretical constraints to the existence of 
a-oscillations. However, as obtained from our analysis, the a-mode fraction of the 
total EWSR diminishes with an increasing mass number. Thus, we do not expect 
that superheavy nuclei will have any relevant a-dipolar response. 
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9 


Conclusions 


We reviewed the phenomenological model of an a-particle moving in a pocket-like 
quasimolecular potential, explaining the gross features of the decay data between 
ground as well as to excited states. We also extended the Brink model to describe 
the a-cluster dynamics in a similar way to giant resonances, thus evidencing the 
possibility to detect a-like collective resonances. 
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